Quantum phase transition in the two-band Hubbard model 
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The interaction between itinerant and Mott localized electronic states in strongly correlated ma- 
terials is studied within dynamical mean field theory in combination with the numerical renormal- 
ization group method. A novel nonmagnetic zero temperature quantum phase transition is found 
in the bad-metallic orbital-selective Mott phase of the two-band Hubbard model, for values of the 
Hund's exchange which are relevant to typical transition metal oxides. 
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The nature of the metal insulator transition in 
multi-component systems, with competing energy scales 
associated with differing bandwidths and Coulomb 
and exchange interactions, is currently of great inter- 
est For instance, in transition metal oxides such 
as Ca2-i:Sr:i.Ru04 or manganites (e.g. Lai-i^Sr^^MnOs), 
the Coulomb interaction between weakly and strongly 
correlated subbands can give rise to complex phase 
changes as a function of temperature, pressure, or impu- 
rity concentration. An interesting aspect of such systems 
is the coexistence of itinerant wide band electrons with 
partially or completely localized narrow band electrons. 
This gives rise to bad- metallic behavior as observed, for 
example, in the resistivity of the paramagnetic phases 
of VO2 and Ca2-xSr:rRu04 at x = 0.2 0. Even 
in standard Mott insulators such as V2O3 and LaTiOa, 
the presence of inequivalent orbitals leads to strong or- 
bital polarization where in the vicinity of the transition 
nearly localized electrons coexist with weakly itinerant 
bands [3, H, 0|. In the case of cuprates, one of the 
most fascinating aspects is the coexistence of strongly 
and weakly correlated regions of the Brillouin zone, giv- 
ing rise to so-called hot spots and cold spots, with break- 
down of Fermi-liquid behavior in the former and Fermi- 
liquid properties in the latter 0, 0, [^. Since the mo- 
mentum variation of the self-energy is associated with 
spatial fluctuations, an intriguing analogy exists between 
multi-site interactions within a single band and single- 
site interactions among inequivalent orbitals. 

A simple model which captures some of the essen- 
tial physics occurring in the systems mentioned above 
is the two-band Hubbard model consisting of narrow and 
wide subbands, coupled via local Coulomb energy U and 
Hund's exchange J [l^, 11 1. As a result of the various 



energy scales contained in this model, its phase diagram 
turns out to be remarkably rich, as shown in detailed 
studies by many authors [l2, [H, Q [H, El, [13, H El, 
20I . 2H . These studies revealed that differing bandwidths 



and a finite Hund's exchange stabilize an orbital-selective 
Mott phase, in which the wide band itinerant electrons 
coexist with localized spins arising from Mott localiza- 
tion of the narrow band electrons. Further increase of 
the Coulomb energy U leads to a Mott transition in the 



wide band whose nature has remained unresolved for the 
important case of anisotropic Hund's exchange of rele- 
vance to many transition metal compounds. 

The aim of the present work is to identify the char- 
acter of the wide band transition of the two-band Hub- 
bard model at zero temperature for anisotropic Hund's 
exchange. We use DMFT [2II in combination with the 



numerical renormalization group (NRG) method [23|, [24 



25 



26l | in order to achieve a sufficiently accurate descrip- 
tion of low-frequency properties. We find that anisotropic 
Hund's exchange yields a continuous quantum phase 
transition as long as J/U > jc, where jc « 0.1. This 
range of J/U values is typical of transition metal oxides 
(e.g., for manganites). In this region, the spectral density 
A{iL!) vanishes at the Fermi level as \Ljf , with S = 1/3. 
This behavior coincides with the one for the metal insu- 
lator transition in the exactly solvable limit U ~ 0, with 
variable J. In striking contrast, the region < J/U < jc 
exhibits a first-order Mott transition, which differs from 
the continuous zero temperature Mott transition for the 
one-band model [2^ . 

We consider an effective low-energy model specific for 
the orbital-selective Mott phase [2lj which allows study- 
ing the transition of the wide band in detail. In this 
model, the electrons in the narrow band of the original 
two-band Hubbard model have undergone Mott localiza- 
tion so their low-energy degrees of freedom are repre- 
sented by localized spins. These couple to the itinerant 
wide band electrons via a Hund's exchange. The result- 
ing model is the ferromagnetic Kondo lattice model with 
interactions in the band, 

Hm = -Y^Ujcl^Cj^ + uY^nii;nii 

ija i 

- Y.[2J Stsl + 2J'{Sts^ +3^3+)] (1) 



where sf and Sj are the z-component and rais- 
ing/lowering operators, respectively, for electrons in the 
wide band at site i and S'f and Sf are the corresponding 
operators for the localized spin at site i. The anisotropy 
of the Hund's exchange is measured by J'/J < 1 with 
J' = J corresponding to isotropic and ,/' = to Ising 
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FIG. 1: (Color online) Effective magnetic moment (s?) per 
site of the conduction electrons for several values of J/f/, 
J' = as a function of U. Solid (dashed) lines: increasing (de- 
creasing) values of U. For < J/U < 0.1 the Mott transition 
is of first order with hysteresis in a region Ud < U < Uc2 (the 
lines are guides to the eye, a finer grid of U points, as used for 
J/U = 0.031, reveals vertical jumps at Uci.2)- For J/U > 0.1 
the transition becomes a continuous quantum phase transi- 
tion. Inset: T — phase diagram and dependence of Uci,2 on 
J/U. OSMP denotes the orbital-selective Mott phase. X de- 
notes the critical point J/U = jc ~ 0.1, U = Uc ~ 1.6 where 
the Mott transition changes from first to second order. 



Hund's exchange. The effective impurity model corre- 
sponding to ll]) which we solve within DMFT-NRG is 
the ferromagnetic Kondo model with a local Coulomb 
repulsion. As in previous works we used a Bethe-lattice 
density of states with half-bandwidth Z? = 1, which 
we henceforth take as our unit of energy. In the NRG 
procedure 800 — 1200 states per energy shell were re- 
tained. The logarithmic discretization parameter for the 
conduction band was A = 1.5 and TV = 90 shells were 
used, corresponding to a minimum energy resolution of 
^-(Af-i)/2 ^YiQ calculations have been carried 

out for fixed values of J/U and increasing values of U /D, 
as we consider the situation where correlations are tuned 
by reducing the bandwidth D whereas the ratio J/U is 
material dependent, being, for example, larger than unity 
for manganites and smaller than unity for most transition 
metal oxides. As we are interested in the metal insulator 
transition we consider a half-filled band and paramag- 
netic solutions only. 

Note that the model ((T]) can also be viewed as a single- 
band Hubbard model in the presence of spin-disorder 
scattering of strength given by the Hund's exchange. For 
[/ = and J > the model reduces to the double 
exchange model used to describe manganites which is 
known to exhibit non- Fermi li quid behavior on restricting 
to paramagnetic groundstates [27[ . This non- Fermi liquid 
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or bad-metallic behavior persists to J7 > 

In the present work we focus on the Ising limit, i.e., 
J' — 0. While this might appear as a severe restriction, 
this limit is actually relevant for all cases with J' < J, 



and therefore should be important in materials where as 
a result of non-cubic octahedral distortions Hund's ex- 
change will not be fully isotropic. This can be under- 
stood from the well known renormalization group flow 
of the ferromagnetic Kondo impurity model [29| and ex- 
plicit numerical calculations including the local Coulomb 
interaction U [s^. For a metallic band, these predict 
that for anisotropic Hund's exchange, J' < J, the spin- 
flip part of the exchange, J', renormalizes to zero, but 
the Ising part, J, renormalizes to a finite value J* > 0. 
Similarly, for an insulating band, one finds that the Ising 
coupling renormalizes to a finite value and the system 
is gapped so J' < J becomes irrelevant for low energy 
properties. From these arguments it follows that, for the 
purposes of determining the nature of the Mott transi- 
tion in the wide band for anisotropic Hund's exchange, it 
suffices to consider just Ising anisotropy J' = 0. Results 
for general J' will be given elsewhere |30| ■ A finite Ising 
exchange for metallic and insulating cases also implies 
that the entropy of the lattice model ([T]) is ln2 per site 
for both metallic and insulating states. Consequently, 
the metal insulator transition at finite temperature need 
not be first order, and in particular a continuous quan- 
tum phase transition can arise at T = which becomes 
a crossover at finite temperature. We also note that the 
renormalization of J' to J* = in the metallic case im- 
plies that the spin-disorder in the lattice model ([1]) per- 
sists to T = 0, giving rise to above mentioned bad-metal 
behavior with Im[S(w = 0)] < [HI, [HI, [111 , where S is 
the wide band conduction electron self-energy. 

To provide an overview of the different transition 
regimes of the model defined in ^ we show first in Fig.[T] 
the behavior of the effective conduction electron magnetic 
moment, (s^), as a function of increasing/decreasing U 
for several fixed values of J/U. For J/U = we recover 
the Mott transition in the single-band Hubbard model 
with a hysteresis region and a continuous metal insula- 
tor transition at U = Uc2 ~ 3. The hysteresis behavior 
persists for small but finite exchange although it rapidly 
diminishes with increasing J/U (see inset to Fig. [T|). 
Within the numerical accuracy of our calculations, the 
hysteresis behavior vanishes for J/U > jc ~ 0.1, indi- 
cating that the Mott transition in this region is a con- 
tinuous quantum phase transition. The critical Coulomb 
interaction for the Mott transition is seen to decrease 
monotonically with increasing J/U (inset to Fig. [T]). 
The precise location and nature of the critical point at 
J/U = jc ~ 0.1, U ~ Uc ~ 1.6 where first and second 
order lines meet requires a detailed study and is out- 
side the scope of the present paper. The phase diagram 
shown in Fig. [1] is consistent with ED-DMFT results at 
low T which also yield first-order behavior for J/U < 0.1, 
U < Uc ~ 1.6 and a continuous transition for J/U > 0.1, 

U >Uc^. 

Let us now look closer at the continuous quantum 
phase transition for J/U > jc- Fig. [H^a) shows the evo- 
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FIG. 2: (Color online) Upper panel: Spectral density of wide 
band for several values of U, for J = 17/4, J' = 0. The Mott 
transition occurs at Uc = 1.183 and is continuous. Inset: 
imaginary part of the self-energy for U < Uc, indicating a 
bad metal. Lower panel: Low-frequency region of spectral 
density A{uj) for J = C//4, J' = at [7^ = 1.183 (red solid 
line; logarithmic plot left inset); the analogous spectrum for 
the exactly solvable model at Jc = 2, (7 = 0, is shown by the 
dashed line. Right inset: ^1(0) as a function of U. 



lution of the spectral density as a function of U for fixed 
J/U =1/4. For J' J, this would be a value typical 
for transition metal oxides which has been used in many 
previous calculations [III, [H, [13, [li, [H, [H, [13, [H, [ij. 
At small U, the density of states at the Fermi level, 
Ep = 0, satisfies the pinning condition, A{0) = 2/7r, 
approximately. Near U ~ I the spectrum begins to ex- 
hibit a pscudogap which gets progressively deeper with 
increasing U. The characteristic scale for this pseudo- 
gap is J. At Uc = 1.183 this pseudogap becomes a 
soft gap before a full insulating gap opens at U > Uc- 
This critical value is in reasonable agreement with ED- 
DMFT results for T = [H and T > [T^. In the 
bad-metallic region below the Mott transition, the elec- 
tronic lifetime t{uj = 0) ~ — l/Im[E(cj = 0)] is finite (see 
inset to Fig. [Ifa)) and vanishes a.t U = Uc- Note that 
for this relatively large value of J/U, well defined upper 
and lower Hubbard bands are only formed once the insu- 
lating gap has opened. In contrast, for smaller values of 
J/U well defined Hubbard bands arc formed already in 




FIG. 3: (Color online) Spectral densities of the wide band 
for several values of U, for, J/U = 0.031, J' = 0. The Mott 
transition is first order and occurs at Uc ~ 2.09. Inset: A{0) 
as a function of U showing hysteresis and discontinuous jump. 
Solid (dashed) lines: increasing (decreasing) values of U- 



thc metallic state, see Fig. [3]below. These T = spectra 
are consistent with those derived within Quantum Monte 
Carlo (QMC) DMFT at finite T [HI, [2^, [sli . Because of 
the presence of the pseudogap the spectra exhibit a char- 
acteristic four-peaked structure, with low-energy features 
limiting the gap and high-energy peaks associated with 
Hubbard bands. 

We have analysed the low-frequency behavior of the 
spectral density at the quantum phase transition in the 
region J/U > jc and find A{ui) ~ \lu\^ with S = 1/3, as 
shown in Fig. ^h)- Moreover, A{0) ~ \U ~ Ud'^^^ for 
U — > U" - Remarkably, this behavior coincides with the 
one obtained for the continuous Mott transition in the ex- 
actly solvable model J/U ^ oo with U ^ on increasing 
J. Solving the DMFT equations for this model we find 
a band-splitting transition at Jc = 2 with A{lu) ^ \lu\^^^ 
and ^(0) ~ I J - Jc|i/2 for J _^ J- ^ This suggests 
that the low-frequency excitations in the whole region 
J/U > jc are governed by the Ising Hund's exchange 
and that the paramagnetic metal insulator transition be- 
longs to the same universality class as the band-splitting 
transition in the model with [7 = on increasing J. 

We now discuss the region J/U < jc where the effec- 
tive two-band model ((T|) yields a first-order Mott tran- 
sition at T = 0. This is evidenced by the discontinuity 
in {sD at both Uci and Uc2 (see Fig. [1)) and in the dis- 
continuous jump in the spectral density at Ep, (see inset 
to Fig. [3|). A comparison of the groundstate energies 
of metallic and insulating solutions in the hysteresis re- 
gions shows that the metallic solution is lower in energy 
than the insulating one, so the first-order transition at 
finite J/U occurs at Uc2 as for the T = metal insulator 
transition in the single-band model. The spectral densi- 
ties for J/U = 0.031 shown in Fig. [3] demonstrate that, 
as for J/U — 1/4, the pinning condition ^(0) = 2/7r 
is increasingly violated with increasing U, indicating the 
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bad-metallic nature of the phase below the transition. 
Because of the smaller value of J the pseudogap is con- 
siderably narrower than for J/U = 1/4. As a result, the 
characteristic four-peaked structure in the bad-metallic 
phase is now better resolved than for the ones presented 
in Fig. [2Ia). Note that, despite the small value of J, the 
spectra shown in Fig. [3] differ markedly from the corre- 
sponding spectra of the single-band Hubbard model ob- 
tained in the limit J = 0. Ising Hund's exchange leads 
to a significant redistribution of spectral weight at low 
frequencies so that even a weak coupling to the local- 
ized spins creates a pseudogap and destroys the zero- 
frequency pinning condition. We find that these effects 
persist also for weak anisotropies J' / J « 0.8 . . . 0.9 [30l |. 

In summary, we have studied the metal insulator tran- 
sition in the two-band Hubbard model with bands of 
different widths and for Ising Hund's exchange, argu- 
ing that the results for the metal insulator transition in 
the wide band apply for all anisotropies J' < J. The 
present NRG-DMFT calculations allow a sufficiently ac- 
curate analysis of the low-frequency properties of the 
wide band in the case of anisotropic Hund's exchange. 
The results reveal a continuous quantum phase transi- 
tion for J/U > 0.1, i.e., in the region relevant for many 
transition metal oxide materials. Close inspection of the 
critical properties of this transition show that they co- 
incide with the ones of the band-splitting transition in 
the model with U ~ Q and variable J. In the region 
J/U < 0.1, the transition becomes first order with hys- 
teresis. This can be interpreted as the stability of first- 
order transitions, in this case that at J = for small 
finite temperatures, to weak perturbations such as the 
weak coupling to localized spins. A sufficiently large per- 
turbation destroys this first-order behavior, giving rise, 
in the present case, to a continuous quantum phase tran- 
sition for J/U > 0.1. The wide band is shown to be 
bad-metallic for all J > 0, with a low- frequency pseudo- 
gap whose width is determined by J and Hubbard bands 
whose positions are determined by U. The resulting four- 
peaked spectral densities differ fundamentally from the 
familiar three-peaked structures of the one-band model to 
which our model reduces in the limit J = [s^ . We have 
also argued that the present results for fully anisotropic 
Hund's exchange (J' = 0) arc relevant to the generic 
case of weakly anisotropic exchange {J' < J) with the 
critical properties of the wide band transition remaining 
unaltered for J' < J. 
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